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We calculate the asymptotic entanglement capacity of the Ising interaction a z ®a z , the anisotropic 
Heisenberg interaction a x <8> o x + o y ® a y , and more generally, any two-qubit Hamiltonian with 
canonical form K = fi x o x ® o x + fj, y a y (g) o y . We also describe an entanglement assisted classical 
communication protocol using the Hamiltonian K with rate equal to the asymptotic entanglement 
capacity. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.67.Hk 



The fundamental resource for quantum information 
processing is an interaction between two quantum sys- 
tems. Any Hamiltonian Hab Ha + Hb that is not 
a sum of local terms couples the systems A and B. To- 
gether with local operations, the coupling can be used to 
generate entanglement [ 110, J| , to transmit classical and 
quantum information , |4|, and more generally, to 

simulate the bipartite dynamics of some other Hamilto- 
nian H' AB and thus to perform arbitrary unitary gates 
on the composite space Hab — Ha ®Hb @, [§[ ||. 

Much experimental effort has been devoted to creating 
entangled states of quantum systems, including those in 
quantum optics, nuclear magnetic resonance, and con- 
densed matter physics [|o|. Determining the ability of 
a system to create entangled states provides a bench- 
mark of the "quantumness" of the system. Furthermore, 
such states could ultimately be put to practical use in 
various quantum information processing tasks, such as 
superdense coding [Tl} | or quantum teleportation 

The theory of optimal entanglement generation can be 
approached in different ways. For example, Ref. [Q con- 
siders single-shot capacities. In the case of two-qubit 
interactions, and assuming that ancillary systems are 
not available, Ref. presents a closed form expression 
for the entanglement capacity and optimal protocols by 
which it can be achieved. In contrast, Ref. (|] considers 
the asymptotic entanglement capacity, allowing the use 
of ancillary systems, and shows that when ancillas are 
allowed, the single-shot and asymptotic capacities are in 
fact the same. However, such capacities can be difficult 
to calculate because the ancillary systems may be arbi- 
trarily large. 

In this paper, we calculate the asymptotic entangle- 
ment capacity of any two-qubit interaction that is lo- 
cally equivalent to [i x a x (S> <J X + \i y a y ®a y , and thus 
present a connection between the results of Refs. |lj and 
0. We consider the use of ancillary systems, and show 
that they do not increase the entanglement capacity of 
these interactions. Thus in these cases, the asymptotic 



capacity discussed in Ref. [2| is in fact given by the ex- 
pression presented in Ref. [ 1 1 . As an application of this 
result, we present an explicit ensemble for entanglement 
assisted classical communication implicitly found in 
Ref. P|, at a rate equal to the entanglement capacity. 
We also give an alternative ensemble achieving the same 
rate. Finally, we conclude by presenting some numerical 
results on the entanglement capacity of general two-qubit 
interactions. 

We begin by reviewing some definitions and known re- 
sults. Let \ip) be a state of the systems A and B. This 
state can always be written using the Schmidt decompo- 
sition |Ol, 



(1) 



where {|</>i)} and {|?7i)} are orthonormal sets of states, 
and \ > with J^. A.; = 1. The entanglement between 
A and B is defined as 



Y^ log A; 



(2) 



(Throughout this paper, the base of log is 2.) 

Reference [|l] considers maximizing the rate of increase 
of entanglement when a pure state is acted on by e~ tHt , 
the evolution according to a time-independent Hamilto- 
nian H (we set K = 1 throughout this paper). We refer 
to this maximal rate as the single-shot entanglement ca- 
pacity. When no ancillas are used, this is given by 



E H :— max 

\ll))(LH A 



lim 



E(e-* Ht \i>)) - E(\i>}) 



t 



(3) 



Here the rate of increasing entanglement is optimized 
over all possible pure initial states of Hab without an- 
cillary systems. In fact, the single-shot capacity may be 
higher if ancillary systems A' and B', not acted on by H, 
are used. For this reason, we may consider the alternative 
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single-shot entanglement capacity 



H ■- 



sup lira 

\i/j)en AA , BB , 



B(e-^|V))-£(|^)) 



(4) 



Note that in Eqs. (||) and (Eh, the limit is the same from 
both sides even though it might be the case that E^*^ ^ 
E^ff in general (and similarly for E^). 

For any two-qubit Hamiltonian H, Ref. Q| shows that 
it is locally equivalent to a canonical form 



Ms > M» > |M*I 



(5) 



In terms of this canonical form, the optimal single-shot 
entanglement capacity of any two-qubit interaction with- 
out ancillas is given by 



E H = "(M* + My) > 



:= 2 max y/ x (1 — x) log 



(6) 

1.9123, (7) 



where the maximum is obtained at xq ~ 0.9168. In addi- 
tion, E^' may be strictly larger than E^*^ when \fi x \ > 

Reference [H considers the asymptotic entanglement 
capacity Eh for an arbitrary Hamiltonian H. Eh is de- 
fined as the maximum average rate at which entangle- 
ment can be produced by using many interacting pairs 
of systems, in parallel or sequentially. These systems 
may be acted on by arbitrary collective local operations 
(attaching or discarding ancillary systems, unitary trans- 
formations, and measurements). Furthermore, classical 
communication between A and B and possibly mixed 
initial states are allowed. Reference [Q proves that the 
asymptotic entanglement capacity in this general setting 
turns out to be just the single-shot capacity in Ref. [Q, 
E H = E$, for all H, so 



Eh = sup 



lim 



(8) 



Note that the definition of the capacity involves a supre- 
mum over both all possible states and all possible interac- 
tion times, but in fact it can be expressed as a supremum 
over states and a limit as t — » 0, with the limit and the 
supremum taken in either order. 

Let \ip) be the optimal input in Eq. (|j) or (||). When 
\ip) is finite dimensional, the entanglement capacity can 
be achieved |i} || by first inefficiently generating some 
EPR pairs, and repeating the following three steps: (i) 



transform nE(\ip)) EPR pairs into 



®7l 



14 if 



(m) 



evolve each according to H for a short time St, and 
(Hi) concentrate the entanglement into n(E(\ip)) + 6tEH) 
EPR pairs 111] . 



In this paper, we show that E^ 



E*£*^ for any two- 



qubit Hamiltonian with canonical form 



K 



tlx 0- x 



fly o~ y 



Mx > My > , (9) 



so 



that all three entanglement capacities are equal: 



E K = E$ = E 



,(i*) 



(10) 



The optimal input is therefore a two-qubit state, and the 
optimal protocol applies. In particular, for these Hamil- 
tonians, which include the Ising interaction a z ® cr z and 
the anisotropic Heisenberg interaction o~ x ® o~ x + o~ y <£> a y , 
entanglement can be optimally generated from a 2-qubit 
initial state \if>) without ancillary systems A'B' . As men- 
tioned above, this result is not generic, since ancillas in- 
crease the amount of entanglement generated by some 
two-qubit interactions, such as the isotropic Heisenberg 
interaction a x ® a x + a y ® a y + a z ® a z H . 

In the following, we will focus on computing the 
asymptotic entanglement capacity of the interaction 



K, 



(11) 



One way to see that this is sufficient to determine the 
asymptotic entanglement capacity of K in Eq. @ is to 
note that K is asymptotically equivalent to 



K' := (fi x + fi y )cr x (g) cr x 



(12) 



and that E t n = for two-qubit Hamiltonians. The 

asymptotic equivalence of K and K' is based on the fol- 
lowing two facts: (i) K' and fast local unitary trans- 
formations on qubits A and B can simulate K |J; con- 
versely, (ii) the Hamiltonian K can be used to simulate 
K' given a catalytic maximally entangled state, without 
consuming the entanglement of A'B' , which subsequently 
can be re-used |J. Therefore, Hamiltonians K and K' 
are asymptotically equivalent resources given local op- 
erations and an asymptotically vanishing amount of en- 
tanglement. In particular, Ek — Ek 1 - This equivalence 
could be generalized to other capacities, but for the spe- 
cific case of entanglement capacity, a simpler proof is 
available. The simulation (i), which does not require 
a catalyst, demonstrates Ek < Ek 1 ■ After computing 
Ek> , we will see that the protocol of Ref. jl| saturates 
this bound, so in fact Ek — Ek 1 with no need for ancillas 
to achieve either capacity. 

We now present the optimization of Eq. (||) for K xx . 
We suppose that in addition to the qubits A and B on 
which K xx acts, c?-dimensional ancillas A' and B' are 
used, where d is arbitrary. We can always write the 
Schmidt-decomposed initial state \ip) as 

= ^2V\\<t>i)AA'^h) B B' (1 3 ) 

i 

= [U® F)(n/A®/ bb 0I $ ) (1 4 ) 
= UVAV T '®7bb'|*) , (15) 
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where U and V are unitary matrices on Haa' and Hbb', 
A is a diagonal matrix with diagonal elements An = A.;, 



m = Ei\i) A A> 



I BB 



i , and we have used the fact that 



(J® M)|$) = {M T <g> I)\<S>) 

for any operator M. Defining p := trBB< 
tanglement capacity of any Hamiltonian H is 



(16) 
the cn- 



Eh 



sup tr 



dp d log p 

--lo g p-p— 



SU!>h ' ~^ l0gP 



(17) 



The variation of p can be computed using perturbation 
theory |: 

^ = -itx BBI [H,\^\]=2lmtx B B'{H\'iP)^\). (18) 

Letting R = U \/~AV T and considering H = K xx , we have 
tv BB > (K xx \iP)(iP\) 

= tr BB ' [(X ® X)(R® iBB'^Mitf ® I BB >)] 

= Hbb' [{XRX T ® Ibb>)\$){$\{R ] ® Ibb')] 

= XRX T R f , (19) 

where we have introduced X := a x ® /, with the identity 
operator acting on the ancilla. The first equality follows 
simply from substitution of K xx and \ip) by their expres- 
sions in Eqs. jll| ) and ([l5]); the second uses Eq. Jlfi|); and 
the third employs the fact that for any operators Mi, M2, 

tr Bfl ,[(Afi <8> /BB/)|$X*l(Af2 ® W)] = MiM 2 . (20) 
bincc p we have 



A'.,. 



sup tr 



_C/t^;7 log A 
at 



(21) 



Using Eqs. ( |18j ) and (|19j), and introducing the Hermitian 
operators Xjj = U'XU and Xy = V'XV, we have 



U^U = 2ImX £/ VAX^VA. 
di 



(22) 



Letting J7, V, A attain the supremum in Eq. ( plf ) (up to 
an amount that can be made arbitrarily small), we find 



E Kmx = -2Imtr (x u \/AX£VAlogA 

(XuVAXy* - Xy VAXjj ) VA log A 



i tr 



itx[M{Xu o X v )] 



(23) 



where we have introduced the real, skew-symmetric ma- 
trix 



M, 



■s/XiXj log(Aj/Ai) • 



(24) 



and the symbol o denotes the Hadamard (i.e., element 
wise) product of matrices. In the second line of Eq. 
we have used 



ImtrA = tr(A- A f )/2i 



(25) 



and the fact that A, Xjj, and Xy are Hermitian. The 
last line can be checked by explicitly writing the trace in 
terms of matrix elements. 

From Eq. ( p3| ) we obtain the following upper bound for 
Ek xx (here \A\ denotes the element-wise absolute value, 



E 



A , 



< tr(\M\\Xu o X v \) 

< suptr(}Af}P) 

p 



< 2 max \J x(l — x) log[x/(l — x)\ 

X 

= a w 1.9123, 



(26) 



where P is a permutation operator and x £ [0,1]. The 
first line uses the triangle inequality. The second in- 
equality follows from noticing that \Xu oXy\ is a doubly 
substochastic matrix |l6[ ]. Indeed, for any two complex 
numbers v and w one has that 2\vw\ < \v\ + \w\ 2 , and 
consequently, for any two unitary matrices V and W , 

E iwai ^ E W + i^'i 2 )/ 2 = L 

i i 

which implies that the matrix \V o W\, with entries 
\VijWij\, is doubly substochastic. Therefore a doubly 
stochastic matrix Q exists such that \Xjj o Xy\ij < Qij 
for all i and j fTJ, so that tr({Mj \X v oX v \) < tr(fMfQ). 
But Q is a convex combination of permutation opera- 
tors Pk, Q = ^fcPfe-Pfej which implies that tr \M\Q < 
supp tr(fMfP). Finally, the third inequality in Eq. ( |2^ ) 
follows from noticing that 



\M\ij = v ^A7|log(A,/A J )| 



< (Ai + A?) max — x) logfx/ (1 — as)] 

a; 

= (Ai + AjOa/2, (28) 



and that 



tv(\M\P) < | E( A ' ; + = « E A * = a ' ( 29 ) 

y * 

where we have used the facts that P is a permutation 
matrix and that A^ = 1. Comparison of Eqs. (0) and 
( p6| ) shows that, indeed, P/f^. = Ejj, completing the 
proof. 
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We have shown that ancillary systems are not needed 
when optimizing entanglement generation by any two- 
qubit Hamiltonian with canonical form given by Eq. (|^). 
More specifically, there is a universal optimal two-qubit 
initial state given by M 



^olO)^ ® \1} B - iVT^\l) A ® |0) fl . (30) 



As an application of the above, we discuss how to use 
the Hamiltonian K to enable classical communication 
between Alice and Bob. This has been studied in 0, 
in which the entanglement assisted forward classical ca- 
pacity (the maximum rate for the Hamiltonian H 
to communicate from Alice to Bob when free, unlimited 
shared entanglement is available) is shown to be 




FIG. 1: Numerically optimized entanglement capacity of the 
two-qubit Hamiltonian p x <t x ®(t x + fi v a y ® a y + a z ® a z with 
single qubit ancillas on each side. The vertical axis in the left 
figure is in units of a. 



Ct{H) =sup 

£ 



lim 

t->o 



x(tr AA , e~ mt £) - X (tr AA , £) 



f 



,(31) 



where £ = {pi, \ipi)} is an ensemble of bipartite states, 
e -iHtg anc j tr^A' £ denote the respective transformed 
ensembles {p u e~ lHt \^i)} and {Pi,ti AA < l^iX^il}, and 



x({PhPi}) ■= S i^2,PiPi \ ~ ^PiS{p t 



(32) 



is the Holevo information of the ensemble {pi, pi}, where 
S is the von Neumann entropy. Reference || also de- 
scribes a protocol to achieve the rate in the bracket of 
Eq. ( |3l"| ) for any ensemble £. 

For any two-qubit Hamiltonian H, Ref. ||] constructs 
an ensemble with communication rate Eh, which implies 
C^(H) > Eh- This ensemble, which is not necessarily 
optimal, is defined in terms of an optimal state for entan- 
glement generation. This ensemble £\ can now be made 
more explicit for Hamiltonian K in light of our findings: 



Pi ■= l^i) := \^o"|0) A ®|l) B + iy/l-x o \l) A ®\0) B , 



P2 ■■= -, \*h) := Vxq\0)a®\°)b - Wl-x \l) A ®\l) B , 

where xq is defined after Eq. (0). For ensemble £\ we 
find 

x(ix A £i) = S(I/2) - S(ix A |V>i)(V>i|) = 1 - £(lVw» 
X {tx A {e- mK £ x )) = 1 - [£(|Vw» - StE K ] (33) 

and therefore the net rate at which classical bits are 
transmitted is indeed Ax/ St = Ek- 

Next we present an alternative ensemble £2 of product 
states with the same communication rate: 



P2 = \ '^ >:= ^(-ly. 



W 1 ~ x J B 



Here, we use K to simulate K 1 |l, under which the en- 
semble evolves. For ensemble £2, S(tr A |V'*XV , *I) — 0, so 

x(tr A £ 2 ) = #2(2:0) 



£2)) = H 2 (x - 28t^Jx^{T^Xo) 



H 2 {x Q ) + E K 5t 



(34) 



(where H2 is the binary entropy) . Thus the communica- 
tion rate is again Ax/ St = Ek- 

The main difference between these two ensembles is 
that the states in ensemble £\ are entangled whereas the 
states in ensemble £2 are not. In the first case the in- 
teraction K is used to decrease the degree of entangle- 
ment between Alice and Bob or, equivalently, to make the 
states of Bob's ensemble tr^ £\ less mixed and thus more 
distinguishable. The same increase of distinguishability 
for the pure states of Bob's ensemble tr^ £ 2 is achieved 
by conditionally rotating them with K, in a way that 
they become more orthogonal to each other. We note, 
in addition, that ensembles £\ and £2 can be prepared 
using different remote state preparation techniques p"7| . 

In conclusion, we have computed the asymptotic en- 
tanglement capacities of all two-qubit Hamiltonians that 
are locally equivalent to \i x a x ® a x + p y a y <8> a y by 
showing that this capacity can be achieved without the 
use of ancillas. However, as discussed above, ancillas 
are necessary to achieve the capacity in general. Al- 
though we do not have a closed form expression for the 
capacity of an arbitrary two-qubit Hamiltonian, we can 
present partial results in this direction. The numeri- 
cally optimized entanglement capacity of a general two- 
qubit Hamiltonian is shown in Fig. |l| Numerically, we 
find that the optimum can be achieved with single-qubit 
ancillas on both sides. For Hamiltonians of the form 
K^ y = p xy {(T x <£> a x + a y ® a y ) + a z <g> a z , we conjec- 
ture that the entanglement capacity is given by 



E 



K„ 



:{V-PiP2log(pi/p 2 ) [sin^ + fi xy sm(ip - £)] 
+ log(p 2 /P4) [sin<^ + p xy sm(6 - £)] 

+ log(pi/p 4 ) P-xy sin£} (35) 
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where the maximum is taken over p\ > 0, p 2 > 0, [4] 
p 4 = 1 - pi - 2p 2 > 0, and 8,ip,^ S [0, 2n). This ex- 
pre ssion was found by investigating the structure of the \SL 
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nu merical optimum, and it agrees well with the numor 

ical results. It does not seem possible to simplify this 
expression further, which suggests that in general, ca- 
pacities may not have simple closed form expressions, but 
can only be expressed as maximizations of multivariable 
transcendental functions. Nevertheless, it would be use- 
ful to show that this maximization can be taken over a 
finite number of parameters by proving an upper bound 
on the dimension of the ancillas. 

We thank Aram Harrow, Patrick Hayden, and John 
Smolin for interesting discussions. We also thank Michael 
Nielsen for comments on the manuscript. AMC received 
support from the Fannie and John Hertz Foundation. 
DWL was supported in part by the NSA under ARO 
Grant No. DAAG55-98-C-0041. FV thanks John Preskill 
and the Caltech IQI for their hospitality. GV is sup- 
ported by the US National Science Foundation under 
Grant No. EIA-0086038. This work was supported in 
part by the Cambridge-MIT Foundation, by the De- 
partment of Energy under cooperative research agree- 
ment DE-FC02-94ER40818, and by the National Secu- 
rity Agency and Advanced Research and Development 
Activity under Army Research Office contract D A AD 19- 
01-1-0656. 



[1] W. Dtir, G. Vidal, J. I. Cirac, N. Linden, and S. Popescu, 

Phys. Rev. Lett. 87, 137901 (2001). 
[2] C. H. Bennett. A. Harrow . D. W. Leung, and J. A. 

Smolin, quant-ph/0205057 
[3] P. Zanardi, C. Zalka, and L. Faoro, Phys. Rev. A 62, 

030301(R) (2000); B. Kraus and J. I. Cirac, Phys. Rev. 

A 63, 62309 (2001); M. S Leifer, L. Henderson, and N. 

Linden, quant-ph/0205055 . 



[6] 
[7] 



020510C , to appear in Phys. Rev. A 



[9] 

[10] 
[11] 

[12] 



[13] 
[14] 
[15] 

[16] 
[17] 



D. W. Berry and B. C. Sanders, |quant-ph/0205181 . 
J. L. Dodd, M. A. Nielsen, M. J. Bremner, and R. T. 
Thew, Phys. Rev. A 65, 040301(R) (2002); P. Wocjan, 
D. Janzing, and Th. Beth, Quantum Information and 
Computation 2, 117 (2002); N. Khaneja, R. Brockett, 
and S. J. Glaser, Phys. Rev. A 63, 032308 (2001); G. 
Vidal and J. Cirac, Phys. Rev. A 66, 022315 (2002); P. 
Wocjan, M. Roetteler, D. Janzing, and Th. Beth, Quan- 
tum Information and Computation 2, 133 (2002), Phys. 
Rev. A 65, 042309 (2002); M. A. Nielsen, M. J. Bremner, 
J. L. Dodd, A. M. Childs, and G M. Dawson, Phy s. Rev. 
A 66, 022317 (2002); H. Chen, |quant-ph/0109115| ; G. Vi- 
dal, K. Hammerer, and J. I. Cirac, Phys. Rev. Lett. 88, 
237902 (2002); LI. Masanes, G. Vidal, and J. I. Latorre, 
Quantum Information and Computation 2, 285 (2002). 
C. H. Bennett, J. I. Cirac, M. S. Leifer, D. W. Leung, 
N. Linden, S. Popescu, and G. Vidal, Phys. Rev. A 66, 
012305 (2002). 

G. Vidal and J. I. Cirac, Phys. Rev. Lett. 88, 167903 
(2002). 

Special issue, Fortschr. Phys. 48 No. 9-11 (2000). 

C. H. Bennett and S. J. Wiesner, Phys. Rev. Lett. 69, 

2881 (1992). 

C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. 
Peres, and W. K. Wootters, Phys. Rev. Lett. 70, 1895 
(1993). 

A. Peres, Quantum Theory: Concepts and Methods, 
Kluwer (Dordrecht, 1995). 

C. H. Bennett, H. J. Bernstein, S. Popescu, and B. Schu- 
macher, Phys. Rev. A 53, 2046 (1996). 

H. -K. Lo and S. Popescu, Phys. Rev. Lett. 83, 1459 
(1999); H.-K. Lo, Phys. Rev. A 62, 01231 3 (2000); P. 
Hayden and A. Winter, quant- ph/0204092 ; A. Harrow 



and H.-K. Lo, |quant-ph/0204096 
R. Bhatia, Matrix Analysts, Springer- Verlag (New York, 
1997). 

Reference 0] cites a method due to P. Shor that can be 
used to prepare £ i , and £2 can be prepared by a technique 
due to P. Shor and A. Winter (in preparation). 



